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This paper presents a robust coordinated control of a flywheel energy storage 
system FESS and a doubly-fed induction generator DFIG based wind turbine 
used to smooth the wind-induced output power fluctuations. The overall 
system control combines field and voltage oriented control schemes and 
nonlinear backstepping approach applied first, to the DFIG side converter 
DSC to regulate the stator active power in order to extract the maximum 
power from the wind velocity, and reactive power in order to maintain, at the 
point of common coupling PCC, a unity power factor. Then, to the electrical 
network side converter NSC to regulate DC bus voltage. Finally, to the 
flywheel side converter FSC allowing the storage unit to serve as a buffer 
that stores energy in the case of excess power and retrieves it back into the 
output in the case of power deficiency. Numerical simulations using 
Matlab/Simulink software demonstrate the validity of the proposed control 
strategies in terms of dynamic response, improvement of output power 
quality and robustness against model parametric variation. 
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1. INTRODUCTION 

Environmental awareness, unstoppable rise of energy consumption along with late developments in 
the field of electrical and power electronics engineering are the main reasons why the world is now heading 
toward green, abundant and constantly replenished energy sources. Among a variety of these renewable 
sources, wind power has arise as one of the most economical and viable options for clean electricity 
generation. As this type of power generation continues flourishing worldwide, and as we keep on increasing 
its penetration level, the power system operator will inevitably confront some serious stability challenges. 
Being that the stochastic and fluctuating nature of the wind profile that is amplified in the output power 
(reaching about hundreds of kW in seconds) induces, in case of high wind power penetration or an isolated 
power system, frequency fluctuations as well as voltages fast (flickers) variations. That could considerably 
affect the electrical network operation and stability [1]. Therefore, the WECS must quickly take part in the 
ancillary services rather than only being considered as passive charges from the power system viewpoint [2]. 

Various configurations have been reccomended in order to reduce output power variation. A 
practical solution consists on allocating a short-term energy storage device such as superconducting magnetic 
energy storage (SMES) [3], supercapacitors [4] or flywheel energy storage system (FESS) [5], as an 
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additional reserve that stores the energy in the case of excess power and retrieves it in the case of deficiency. 
Reference [6], presents a comparative study and a thorough full overview of the technical and economic 
performance on the available energy storage technologies. For instance, with their high efficiency, both high 
power and energy densities, and longer cycling life, the FESSs offer a very competitive option over the 
alternatives [7]. 

Among many other relevant works, in [2] a FESS is associated with a permanent magnet 
synchronous machine (PMSM) based wind power plant in order to smooth and improve the output power 
quality. Classical PI control laws along with a fuzzy logic supervisor are used for the power flow 
management between the FESS and the electrical grid while considering the storage device state of charge 
(SoC). Few have been the studies that suggested a profitable application of FESS to a DFIGURE The author 
in [8] has presented a solution that consists on connecting the FESS to the DC bus through an AC/DC 
converter in order to study the effects of reactive and active power quality improvement on load voltages 
using PI controllers. However, traditional control schemes do not always provide pleasing results, since they 
only guarantee the appropriate tracking response at a certain operating point for which the control law was 
already designed. Whereas, when a parametric variation is registered, which is usually the case due especially 
to modelling uncertainties and some physical phenomena, the performance of Pi-based control scheme risk to 
be insurable and gravely deteriorated. 

Accordingly, in attempt to meet higher performances, this paper presents a coordinated robust 
control of a FESS associated to a DFIG-based wind turbine using a systematic and recursive methodology 
used for nonlinear systems control known as backstepping (BS) [9]. Owing to this technique, the feedback 
control law that guarantees the system overall behavior and stability can be efficiently elaborated. Figure 1 
shows the basic configuration of the system, where the stator windings of the DFIG are straightforwardly 
linked to the grid, while its rotor windings are linked via a back-to-back converter. The DFIG side converter 
DSC is used to control the generated reactive and active powers following a maximum power point tracking 
(MPPT) algorithm. The electrical network side converter NSC is used to keep constant the DC bus voltage 
regardless of the direction and magnitude of the power. Finally, the flywheel side converter FSC is controlled 
in order to smooth the wind-induced fluctuations and to improve the quality of the output active power. 
Simulation results performed in Matlab/Simulink verify the performance of the proposed control schemes, in 
terms of output power smoothing, dynamic response and robustness against parametric variation. 



The rest of the paper is arranged as follows: Section 2 presents the modelling of the wind turbine, 
the DFIG and the FESS. The proposed control strategies of the MSC, the GSC, and the FSC are discussed in 
details in Section 2. Simulation and interpretation of the results in Section 3. Finally, a conclusion in 
Section 4. 
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2. DYNAMIC MODELLING OF THE STUDIED SYSTEM 
2.1. Wind turbine model 

The aerodynamic power harvested by the wind turbine P h can be expressed as function of the power 
coefficient C p as follows [10-12]: 

P h =C p (X,p)-P wind =0.5-C p (X,p) P Av 3 , X = R ' Q t/ ( 1 ) 


Where v is the wind speed, p is the air density and A represents the swept area, R is the blades 
radius and Q r is the turbine rotational speed. The coefficient C p represents the aerodynamic efficiency of the 
turbine, which approached in this work by the expression bellow [9]: 
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Where, p the angle of orientation of the blades and X is the tip speed ratio. Moreover, the aerodynamic torque 
C h is obtained from P h by [9]: 
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( 3 ) 


Neglecting the elasticity, friction and energy losses. The gearbox that is used to adapt the rotational speed 
Q r to that of the generator Q m can be modelled using G a gain ratio as follows: 

Q mec =G-n x ; C h = G-C g ( 4 ) 


Finally, the model is completed by the fundamental equation of dynamic [13]: 


dfl 


dt 
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( 5 ) 


2.2. DFIG Model 

The DFIG model is established using the dynamic equations of the stator and rotor currents, 
voltages, and fluxes represented in a two-phase dq reference frame [9, 14]: 
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Where, i s> i r , v s , v r , & s and & r being respectively the stator and rotor currents, voltages, and 
generated fluxes. co r and oo s represent the angular speed of the rotating field associated respectively to the 
rotor and to the stator (rad/s). Similarly, R s , L s , R r , and L r are stator resistance and incuctance and rotor 
resistance and inductance respectively. Whereas, M represents the mutual inductance. The reactive and active 
power as well as electromagnetic torque can be given by, with p being the pair of poles: 
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2.3. FESS modeling and sizing 

FESS operates as an energy reserve that stores the energy in a kinetic form in a high-speed spinning 
mass coupled to the shaft of an electrical machine that provides an electromechanical interface between the 
flywheel and the system [1]. During the charging process, the machine works as a motor to speed up the 
rotational motion of the flywheel and therefore increasing the stored energy, then, the flywheel remains in the 
standby mode. Whereas, on-demand, the machine serves as a generator and decelerate the flywheel retrieving 
the stored energy into the system. The amount of energy stored in a flywheel Efw depends on the mass and 
the form (represented by flywheel inertia J fw ) and the rotational speed Q fw [6]: 

Eftv =0.5-4^, (8) 


However, for safety reasons, the accessible kinetic energy Ef Wl is limited by Qf W -max an d ^fw-min which 
are respectively the maximum and minimum speed of the flywheel [15], giving that: 

Efwl “ ' Jfw (^fw-max “ ^fw-min ) j (9) 


The flywheel inertia is determined using the expression below, after selecting the storage system maximum 
power capacity stored during a certain At, which is identical to the electrical motor rated power P r _ /M . Which 
in this application, considered lower than 20% the DFIG based wind turbine rated power: 
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In this work, the flywheel exchanges the energy with the system by means of an induction machine 
(IM). It is essential for the IM to operate in its field-weakening region, at a rotational speed ranging between 
1500 and 3000rpm, where its rated power is always available [1]. Therefore, the regulation of the flux cp rd is 
conducted in accordance with the following control law: 
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( 11 ) 


Where, cp r - r ated an( J Vs-rated are respectively the rotor and the stator rated fluxes, Q m and Q r being the 
measured and the rated rotational speed of the IM, and V s is the stator rated voltage. 


3. PROPOSED CONTROL STRATEGIES 
3.1. MSC control 

In order to decouple the control of stator reactive and active powers, a field oriented control scheme 
is adopted, aligning the d-axis with the stator flux (& sd = & s ; & sd = 0). Moreover, by neglecting the stator 
resistance, which is a rational approximation for high power generators [9], we determine the state model of 
the DFIG as follows: 
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The backstepping (BS) basically uses the so-called (virtual control variables (VCVs)) to 
systematically dissacociate a nonlinear complex control scheme into simpler design problem [9]. The control 
design is divided into multiple steps. In every step, we virtually deal with a single input single output 
subsystem, so that each step provides a reference for the next one. The control law that guarantees the overall 
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system behaviour and stability can be built up in a constructive manner using Lyapunov functions. As it is 
shown in Figure 2 The block diagram of DSC control involves an open loop MPPT [16] used to establish the 
optimal active power reference in order to maximize the converted power by adjusting the wind turbine 
mechanical to its optimal value, and two serial control loops. The outer loop provides the VCVs (the 
references of currents) via the reactive and active power control. While the inner loop provides the actual 
control variables (rotor voltages) [17, 18]. 


GRID 



a) Currents references computation 

First, we define the errors variables, with P s _ r and Q s _ r being respectively the references of the stator 
reactive and active power. Then we select the chosen Lyapunov functions (LFs) such that: 


e l P s-r P s^ e l Ps r Ps -—q 2 V --Q 2 

5 M. — 2 e l V 2. “ 2 62 

e 2 = Qs-r - Qs -> e 2 = Q s _ r - Q s 

The derivative of Lyapunov function can be written: 

Vi. =e i-e 1 =e 1 [^Ps-r-(a(v rq +R r I rq +TiI rd -x)j 
V 2 . =e 2 -e 2 =e 2 Q s _ r -^i(V rd +R r I rd -r|I rd ) 


(13) 


(14) 


According to the BS theory, V x and V 2 must be strictly negative to ensure the stability of this first step. Thus, 
Yi = -b sl e^ and V 2 = — b s2 e 2 . The currents references are obtained using (14): 


= ^(^ s - r+bsiei )-i +1 ^ -4 

trd-r = ^-(Vr + bs 2 e 2 j ~ - Tll rd ) 


(15) 


b) Rotor voltages computation 

We define the errors variables, noting that the VCVs determined by (15) are the desired variable in 
this step, and we select the LFs for this step such that: 


e 3 “ f rq-r - bq e 3 ~ Lq-r - Irq 
e 4 = bd-r — bd ^ 64 — Ird-r — Ird 


l V 3 .=bef+e 3 2 ) V 4 .={(e 2 +e 2 ) 


(16) 
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In order to insure a stable tracking performance the derivative of the LFs must be negative. Finally, owing to 
(16) the control voltages are given by: 


-a) 


Vi=-b sl ef-b s3 e^0 
V 4 = —b s2 e^—b s4 e^0 


;-b) 


V rq-r = j- b s3 e 3 -bq-r-^e, [ -{R r I rq + r|I rd + x} 


V rd-r = | — b s4 e 4 — Ir d -r — “ { R drd “ ^rd } 


(17) 


Noting that b s i, b S 2 , b S 3 and b S 4 are positive constants called the BS coefficients, used to ensure a fast dynamic 
response and the stability of the overall system. 

3.2. GSC control 

The GSC is connected to the grid through a filter (resistor-inductor circuit R f , L f ). Figure 3 shows 
the block diagram of NSC control that involves two consecutive control loops. The outer loop provides the 
VC Vs (the references of the filter output currents (7/cZ- r and lf q - r ) via the reactive and active power control. 
While the inner loop provides the actual control variables (the filter output voltages t^ d _ r and Vf q _ r )). This 
system can be represented in a dq reference frame by the following equation [19]: 

v td.= RJtd + L t ~~ _ ^sMtq + v sd 

* ( 18 ) 

V tq . = R,I tq +L t ^-(o s L t I td +V sq 

A vector control is used to achieve an independent control of the active power and the reactive 
power and exchanged between the NSC and the grid. By aligning the q-axis of the park reference frame with 
network voltage equivalent vector, thus, V sd =0 and V sq . = V s . , the model of the system can be determined 
as follows: 


k = ^(v tq ~R t xI, q -a> s L t I td -V s ) 
hq- = T (%q — P-t X hq — thd — ) 

Z l. = P t. 


Q..=^(v, d -R«xI td -a>sL,I« q ) 

■ !td. = r— (V td - R t x I td - co s L t I tq ) 

Z 2. = Qt. 


(19) 



Figure 3. NSC control 


a) Currents references computation 

We define (e 1? e 2 ) representing the errors variables, we choose the FFs such that: 
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6] — ?t-r — ^ Gl — Pt-r — Pt 
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Likewise, applying the same principle used for the DSC control, the VCVs that guarantee a stable tracking in 
the first control loop can be determined such that: 


Vi = -b^,ef -< 0 
V 4 = -b^e^ -< 0 


;-b) 
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( 21 ) 


b) Rotor voltages computation 

The VCVs determined previously are the desired variable for this step. We define : 

e 3 = Itq-r _ \ e 3 = hq-r “ Lq _ 

£4 = Itd-r — ltd ^ ^4 = Itd-r — ltd 

Similarly, the control law that ensures global stability of the system is written such : 


V 3 = ~( e l 2 +e 3 ) 

V 4 =fH + eJ) 


( 22 ) 
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(23) 


c) DC-link viltage control 

Neglecting the converters losses. With P r , P c and Pst being respectively the rotor, the DC-bus and the 
FESS active power Figure 4 We have: 


P c =P t -P s -P st 

with: Pc= i(i cu ^ (t) ] = CUdc(t) '^dT^ 


According to (24) and (25), we write: 

j =c dU dc (t) P t -P r -P st 
c dt U dc (t) 


(24) 

(25) 


(26) 



Figure 4. DC-link voltage control 
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3.3. FSC control 

The FSC control it is used to manage the storage device operation so as to store the energy when 
available, and to restore it back to the system in the opposite case. The block diagram of FSC control is 
shown in Figure 5 it involves two control loops i.e. the outer loop that provides the VCVs (the references of 
the stator output currents / 5d _ r and l sq - r ) via the rotational speed and direct flux control, and the inner loop 
that provides the actual control variables (the output voltages F scZ _ r and V sq _ r )). 



Figure 5. Block diagram of the FSC control 


The reactive and active power at the point of common coupling (PCC) are given by: 

(27) 

(28) 

Where (Q s ; P s ), (Q r ; P r ) and (Q t ; P t ) being the reactive and active power provided in the stator, the 
rotor and the filter, respectively. Knowing the adjusted power P g _ r desired to be injected into the grid, the 
power reference of the FESS control can be determined using the following expression [6]: 

Pstk=Pg-r-P m (29) 

Therefore, with E fw0 being the energy initially stored in the flywheel. We determine the total power stored in 
the flywheel using: 

Efiv-r =E fw o+ Jp st k dt (3°) 

Finally, the speed reference can be deduced using the expression (11): 


p g =p s +p, 

Q g = Q S +P t 

p =p +p 

m c 1 r* 


We define: 


:Q S +Pr 


Q 


fw-r 


| 2-E fa _ r 
V J fw 


(31) 


In order to achieve a separate control the rotational speed and the flux component of the IM, we adopted a 
field oriented control scheme, by aligning the d-axis of the park reference frame with rotor flux, such that 
cPrd — <Pr> <Prq = 0- Based on the IM model given in [20] we determine the stat model of the IM as follows: 
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Q = ±^ (Prlsq-— n ~ — 
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isq =-r|I sq -ppai sd -pQI sd -aM^ik-^ 3 - 

<Pr ctL s 


y, =Q 


<Pr =-a<p rd +aM I sd 

Isd =-riI s d+ a P(Prd+P nI sq- 

y 2 = <Pr 


( 32 ) 

<Prd aL s 


Where, p represents the pole pairs and: 


t M 

v ^ r L s j 


pM M 2 R r +L 2 R r n M R r 

r ; r|= — t2 ; P=-j-i“ ; a= ~T 

L aL c L r crL L L 


a) Currents references computation 

We define (e 1? e 2 ) representing the errors variables and we select the lyapunov function such that: 


e, = Q r - Q —» ei = Q r - Q 
^2 = 9 r-ref "<Pr e 2 = <Pr-ief-<Pr 



(33) 


Likewise, applying the same principle used for the MSC-GSC control, the VC Vs that guarantee the 
asymptotic stability of the first control loop can be determined such that: 
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vi=-b; ie 2 ^o 
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; -b) 
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. \ 
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(34) 


b) Rotor voltages computation 

Again, the VC Vs determined previously are the desired input variable in this step. We define the errors 
(e 3 , e 4 ) and Lyapunov function such that: 

e 3 = Irq-r ~~ ^rq e 3 = Irq-r - Irq 
e 4 = I r4 _ r — I r( j —> e4 = Ird-r — Ird 

The control law that insures the stability of the overall system, can be determined by: 

V sq-r = CTL sfVr + Y<P r e 1 -5 1 +b" 3 e3 

„ ^ 

V sd-r = aL s Ird-r + aMe 2 -8 2 +b s4 e 4 


II 1 

Where: 8j = -r|I - p(iO<p, - pQI sd - aM 5 2 = -r|I sd - a(3<p r - pQI -a.M^ 

<Pr <Pr 


-a) 


V3=-b: i e?-b: 3 e 2 ^0 
V 4 = -b' 2 e 2 -b' 4 e 4 -< 0 


;-b) 




V, =-|e?+e 2 ) 


V 4 =i(e 2 + e 


(35) 


4. SIMULATION RESULTS 

Simulation studies are performed with the system characteristics presented in Table 1 in Appendix. 
All simulation are run with the same wind profile that varies around an average speed of 8.9 m/s for a 50 sec 
duration, as shown in Figure 6. 
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When using a MPPT strategy, the active power delivered by the wind turbine is corresponding to the 
cube of the wind speed. Meaning that a small variation in the wind causes a significant variation in the output 
power. Figure 7 shows the DFIG stator active power and and Figure 8 shows its reactive power curves 
respectively. In this work, the latter is maintained equal to zero in order keep a unitary power factor at PC. In 
brief, the results demonstrate the validness of the MSC proposed control scheme, being that the outputs 
encounter flawlessly their references and all tracking errors converge to zero asymptotically. Moreover, the 
reference of the FESS power P stk-r determined using Eq.42, which corresponds to the difference between 
P g _ r which is the net power intended to be delivered to the grid and P m the total generated power, as it is 
shown in Figure 9. In Figure 10 and Figure 11, the results show the effectiveness of the FSC control scheme, 
given that the FESS output power P stk and the instantaneous rotational speed Q present an excellent tracking 
response. The storage device is initially charged with E fw0 and rotate at an initial speed of 200 m/s, it 
accelerate to store energy when the output power P stk is positive, and decelerate to generate it in the opposite 
case. Figure 12 and Figure 13 pinpoints the efficiency the GSC control strategy. Figure 12 illustrate the DC 
bus regulation, which is used to maintain a constant voltage, and Figure 13 shows that the reactive power Q f 
exchanged between the GSC and the electrical network converges to its reference that is likewise chosen to 
be zero. 

The total active power injected to the power system P g with and without including the FESS support 
is shown in Figure 14. As expected, the power injected without the support of the storage device is highly 
fluctuating due to random nature of the wind. On the other hand, the quality of the output power generated by 
the arrangement (DFIG plus FESS) at the PCC is comparatively smoother, and P g remains at the designated 
reference P g _ r which chosen to be equal to 775kW in this case of study. However, a low pass filter can be 
used to determine P g _ r for long-term FESS operation. 

In order evaluate the improvement achieved by applying the proposed control algorithms, it is 
particularly important to analyze the system performance under parametric variation, since usually machines 
parameters are exposed to many inaccuracies due to the adopted identification methodology, measuring 
devices or several natural phenomena. Therefore, to demonstrate the system ability to react against this kind 
of disturbances, we compared the proposed backstepping control approach with conventional PI controllers 
for the MSC and the FSC control. This test, consist of modifying the DFIG’s model parameters such that 
(R' r = 1.5R r ; L' r = 1.2L r ; L' s = 1.2L S ) and the FESS-IM’s model parameters in a way that (R' s ' = 
1.5R S ; R' r = 1.5R r ). Figure 15, Figure 16 and Figure 17 show respectively the DFIG stator active P s and 
reactive power Q s and the FESS output power. Accordingly, the test confirm the robustness of the proposed 
backstepping control strategy compared to PI controllers. Being that the latter have shown high sensitivity 
towards parametric variation, a static error appears in the FESS-IM control, notably at the output power, 
whereas the performance has been completely deteriorated for the DFIG control. On the other hand, BS 
controllers present an exquisite disturbance rejection and all the outputs converge correctly to their 
designated references. 



Time (sec) 

Figure 6: Wind speed profile 
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Figure 12. Reactive power at the output of GSC 
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Figure 7. Active power of the DFIG 


Figure 13. DC-link votage 



Figure 8. Reactive power of the DFIG 


Figure 14. Power injected to the grid 



Figure 9. Operation principle 



Figure 10. Output active power of the FESS 



Figure 11. Mechanical speed of the flywheel 



Figure 15. DFIG Active power (robustness test) 
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Figure 16. DFIG reactive power (robustness test) 



Figure 17. FESS active power (robustness test) 


Int J Pow Elec & Dri Syst, Vol. 10, No. 2, June 2019 : 1110-1122 













































































































Int J Pow Elec & Dri Syst 


ISSN: 2088-8694 


□ 1121 


5. CONCLUSION 

In this paper, an IM based FESS was associated to a DFIG based wind power plant in order to 
address the issue of output power fluctuations. A novel robust control algorithm based on field-oriented 
control schemes and nonlinear backstepping approach was applied, first, for the DFIG control through the 
DSC in order to ensure MPPT operation and a unity power factor at the PCC. Then, to the NSC for DC-bus 
voltage and reactive power (exchanged between the NSC and the power system) regulation. And finally, to 
manage the storage unit operation through the FSC so as to store the energy when available and retrieves it to 
the power system in the opposite case. Simulation results have shown in one hand, the enhancement 
introduced to the quality of the power injected with the support of the storage device. On the other hand, they 
validate the improvement achieved by applying the proposed control algorithms compared to classical 
controllers under ideal operating condition and in a case where models parametric variation is registered. 
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APPENDIX 

Wind Turbine parameters. Rotor blade radius R=35m, mechanical gearbox G=60. Moment of inertia J=4M 
kg/m2. Damping coefficient f=0.017 N.m.s/rd and the density of air p=1.2 kg/m3. 

DFIG parameters: Rated Power 1.5Mw, grid voltage line to line V s =690; f s = 50 Hz, number of pole pair 
p=2, stator & rotor resistances and inductances R s =12mQ, R r =21mQ, L s =13.7mH, L r =13.6mH, respectively. 
Mutual Inductance M=13.5mH. 

FESS (IM) parameters: Rated power 450kw, number of pole pair p=2, stator & rotor resistances and 
inductances R s =51mQ, R r =51mQ, L s =47.1mH, L r =40.7mH, respectively. Mutual Inductance M=40.1mH, 
flywheel inertia J v =729.5 kg/m 2 and the viscous friction coefficient f=0.008 N.m.s/rd. 

DC-link & the filter parameters: R t =2mQ, L t =5mQ, DC link capacitor C=4400pF and voltage Ud C =1750V. 
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